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In the framework of Newtonian cosmology or general relativity it is simple to derive a mass scale 
below which collapsed structures are relatively devoid of baryons. We examine how the inclusion of 
a chameleon scalar field affects this baryonic Jeans mass, bearing in mind both the canonical case of 
a gravitational-strength coupling between the scalar field and matter, as well as the strong coupling 
regime wherein the coupling is very large. We find that baryon effects persist down to smaller scales 
in a chameleon theory than they do in ordinary general relativity, especially in the case of strong 
coupling. Several potentially observable consequences of this are identified. 



Modified theories of gravity often invoke a screening 
mechanism to remain consistent with precision tests 
of gravity on the terrestrial and solar system scale [T}- 
The chameleon scalar field 7, 8 offers one such 
mechanism; in this family of models a scalar field 4> af- 
fects matter through a conformal coupling {(f>) . The 
Einstein-Hilbert action is unchanged from ordinary GR 
but the matter fields are assumed to be coupled to 
a Jordan frame metric g^i, = A^{4>)g^,^. Assuming a 
Minkovifski metric, negligible time evolution of (j) and 
non-relativistic matter of density p, this coupling has 
the effect of introducing a density dependence in the 
effective potential which governs the evolution of the 
field: 

V2</, = aVe///9</)2 (1) 

where 

Veff = V{<j,)+pA{<j,), (2) 

The presence of this term in the potential means that 
for suitably chosen V{(t)) and A{(j)) the ambient matter 
density dictates the mass of the field m? — d^Vef f / dcjP . 
In high density environments like the Earth the mass 
becomes very large relative to low density environ- 
ments like the cosmos where the field could be so light 
as to play the role of dark energy [7J [H] . This ability 
to 'blend in' with its surroundings earns the field its 
name. 

The purpose of this paper is to connect several recent 
findings regarding the chameleon model. The original 
chameleon model assumed A{(j)) = e^"^/*^^', where /3 is 
the coupling constant and is taken to be order unity, 
corresponding to a gravitational-strength coupling be- 
tween matter and the chameleon field. This model 
demonstrates the surprising result that a scalar field 
can couple to matter non-negligibly and yet evade ter- 
restrial constraints on equivalence principle violation 
and deviations from the gravitational inverse-square 
law [3 IH]. However, Mota and Shaw later demon- 



strated fH" that in fact the chameleon could be much 
more strongly coupled (/3 ^ 1) and still evade terres- 
trial detection. 

On the other hand, Brax et al. have studied 
the evolution of dark matter density perturbations 
in chameleon cosmologv|10j. and found that on suf- 
ficiently small scales the perturbations grow with an 
exponent that depends linearly on /3. Even for the 
case of /? ^ 1 small scales demonstrate dramatic en- 
hancement within the chameleon model. Thus we are 
motivated to ask whether the strong coupling regime 
would produce astrophysical signatures in the realm of 
small scale structure. 

The analysis of Brax et al. does not include the ef- 
fect of baryon pressure on the evolution of cold dark 
matter (CDM) perturbations. It is well-known that 
in ordinary GR baryon perturbations follow the CDM 
perturbations on large scales [TT]; due to the acous- 
tic pressure of baryons there is a Jeans scale below 
which collapsed structures have negligible baryon con- 
tent. To illustrate the kinds of effects the chameleon 
might have on small scale structure, here we study the 
same problem in the context of chameleon cosmology 
by generalizing the analysis of Brax et al. to include 
baryons and acoustic pressure. 

First we review the baryonic Jeans calculation within 
general relativity, i.e. the growth of linear perturba- 
tions in CDM and baryonic matter density, from decou- 
pling at redshift Zdec ~ 1100 to the redshift at which 
baryons and radiation are still at the same tempera- 
ture, z = Zf, ^ 150. The presentation here closely 
follows Section 8.3 of Weinberg [T?|. First order den- 
sity perturbations of CDM and baryons {Spo and dps 
respectively), to their otherwise uniform backgrounds 
and jDg obey identical continuity equations 

^^^^ + SHSpD - a^^p^q^SuD = (3) 



+ 3H5pB - a ^pgC^duB = 0, (4) 



2 



at co-moving wave-vector q. However, the Euler equa- 
tion for each type of matter is different, reflecting the 
fact that baryons experience pressure. The baryonic 
pressure is reflected in the speed of sound, denoted 



dps/dpB- Thus, the Euler equations are 

dSuo AnGa 
— — + HdUD = — ^ [opD + OpB) , 
at q-^ 



(5) 



relevant (growing) solution 
1^ = 2/3, $ 



1 



l + a 



It is convenient to write a = q^/qj where 

9 „ 2 M^AT 



(13) 



(14) 



dSuB rrc 47rGa 

—J- \-HduB = — 5— (opD + <5ps) —opB- (6) 

at q'^ apg 

Note that we can eliminate the velocity potentials 5u o 
and 6ub via Eqs ([3| and Q. Writing Sd = Spo/Po 
and 6b = 6pb/pb^ ^^'^ using the fact that both p^ 
and pg fall off as a~^ or t^^ in this era, we can write 
the Euler equations ^ and Q as 



3D 



3*2 



ifSB + (1 - DSd) 



(7) 



and 



•■ 4 • 2iT 2 

^s + 3fe--^fe + ^(/fc + (l-/)M (8) 



where 



3q2i;2t2 



^ ^'q' /=i^ = ik (9) 

2a2 47rGjOMa^' Pm 



and Pjvj = p^, + ^^3. 

Fo the epoch we are considering, roughly 1100 < 
z < 150, the temperature T of baryonic matter was 
the same as radiation and cx T oc a""'^, so a was 
constant. Thus we can seek power law solutions to 
Eqs. ([7| and If we write 



(10) 



Note that the two parameters v and ^ encompass the 
effects we wish to examine when we generalize to the 
chameleon case: the growth of the perturbations and 
the relative strength of the baryon and CDM pertur- 
bations. We assume they are independent of time but 
not wave-vector q. With these assumptions Eqs. ([7| 
and dsl become 



^^ + T = T,{n + {i-f)). 



(11) 



Here mjv is the mass of a hydrogen atom and p = 1.22 
takes in to account that matter was a mixture of hy- 
drogen and helium atoms. i?o is the present value of 
the Hubble constant and T^o is the present temper- 
ature of the cosmic microwave background. In going 
from the expression for q in eq ([9| to the expression 
above we have made use of the formula for the speed 
of sound in a gas v1 = 5kBT/3p,mpf . Plugging in these 
values we find that qj = 384/i (Mpc)~^ corresponding 
to a length scale Aj = ^ir/qj = IQ/h kpc. The Jeans 
mass Mj is then defined as the amount of matter in a 
sphere of radius Aj at the present epoch. Numerically 
Mj « 6 X lO^M©. 

Physically, f > 1 corresponds to baryon enrichment 



and ^ < 1 to baryon depletion. Eq ( 13 1 reveals that ^ 



1 for small wave-vectors q <Si qj and C ~^ large wave- 
vectors q ^ qj. In other words on large scales baryon 
perturbations follow dark matter perturbations but on 
small length scales there is significant baryon depletion. 
Equivalently one can say that collapsed structures of 
mass less than Mj have significant baryon depletion. 

It is straightforward to extend this analysis to in- 
clude a chameleon scalar field. Brax et al. showed that 
on sufficiently small scales the effect of the chameleon 
on the evolution of perturbations is simply to enhance 
G by a factor of (1-1-2/32). q^^ s\xc\y scales Eqs. ([7| and 
([s]) therefore become 

5d + \^D = ^{1 + 2/3') {fSB + (1 - f)5D) (15) 



and 

■i 4 . 

Ob + ^Sb 
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(l + 2/32)5s + T^(/^s + (l-/)<5i,). 

(16) 



Employing the same power law ansatz, we have 



z.2 + ^ = ?(l-H2;32)(/eH-(l-/)), (17) 



+ ^ + f = ^(/+(l-/)/0. (12) 

Eliminating ^ yields a quartic equation for v] Weinberg 
makes the approximation that / = and obtains the 



-^ + ^ + f = ^(l + 2/3')(/ + (l-/)/0- (18) 

We need not resort to the putative smallness of / to 
obtain a solution; the quartic equation for v obtained 
by eliminating ^ can be solved exactly. The solutions 



are 



and 



where 



u+± = -(-1 + v/l + 36n±) 



- v/l + 36n±) 



-1 



"± = ^(a-(l + 2/3^))±7 



(19) 



(20) 



(21) 



7 = V(a - (1 + 2/32))2 + 4a(l + 2^2)(i _ j) (22) 

From this we 



The only non-decaying solution is v- 
det ermine 



2(l + 2/32)(l_/) 



(23) 



a+(l + 2/32)(l--2/)+7' 
We have left / unspecified, which makes it easy to 



check that this reduces to the solution Eq.(13l for 
f — /3 — Q. The baryon fraction to good approxi- 
mation is / = 1/6; for strong enough coupling such 
that 1 + 2(5'^ K 2/3'^, this gives 



-1 1 

T ^ 6 



1 



.36(_i(a_2/32) 

, 1/2 



(24) 



where we have written the terms out explicitly to 
demonstrate the dependence of v on /3. But for present 
purposes we are more interested in the relative strength 
of the baryon perturbation; defining y = q/(3qj, we 
have 



20/6 



y+l + -^iy + 2y-ly' 



(25) 



Note that the above approximation is not so bad 
even for /3 = 1, and the overall behavior of ^ is the 
same, but more specifically for this case we have 

6 = , ^ (26) 

2;2+2 + v/(y^ + 3)2-2y2 

We plot the behavior of ^ for ordinary CDM and 
the two chameleon models in Figure 1. Clearly the 
attenuation of baryon presence with wave number is 
much less rapid in the chameleon case of /3 = 10 than 
for 13 — 1 or ordinary GR. 




FIG. 1: Relative strength of baryon to dark matter pertur- 
bation as a function of wave- vector for ordinary GR/CDM 
(dashed) , and a chameleon model with 13 = 1 (orange) 
and /3 = 10 (green). The /3 = 10 case differs significantly 
from CDM case; baryon perturbations persists at smaller 
structures in a chameleon theory with /3 > 1 relative to 
normal CDM. As discussed in the text, the chameleon only 
affects modes above a characteristic scale qcham', modes 
much greater than this threshold will be influenced by the 
chameleon throughout their evolution and modes much be- 
low this scale will evolve unaffected by the chameleon. For 
modes between these two regimes the chameleon kicks in 
at some point during their evolution. Generally speaking 
Qcham^increases as j3 increases, for the fiducial model in 
Eq.|27|, modes with q > gj/3^''^(70 eV/Mf^^ will undergo 
chameleon-influenced evolution throughout, and those with 
q < qj/3^''^(70 eV/Mf/'^izt/zdecf^^ will not be afi-ected by 
the chameleon for the range of redshifts considered here. 
For the potential in Eq.( |27[ ) it is possible then for a strongly 
coupled chameleon to evade astrophysical constraints by af- 
fecting only unobservably small scales. But for the alterna- 
tive parametrization given in Eq. ( 28 1 , M2 replaces M in the 



threshold wave- vector, and sub-galactic scales and scales of 
order qj would be strongly affected. 



Discussion 

What are the implications of the baryon persistence 
demonstrated in Fig. 1? The index i' represents the 
exponent to which linear perturbations in the matter 
fluid grow with time. Since it acquires a linear depen- 
dence on /3 in chameleon models, things quickly get 
out of control in the strong coupling regime; for ex- 
ample for P = 10, iy{q lOOqj) « 80. Clearly this 
regime needs to be handled more carefully than the 
linear treatment, nonetheless several observations can 
be reliably inferred from the above analysis. 

The chameleon has a characteristic scale kchamSA. 
which it acts; this is a function of redshift as well as the 
particulars of the potential V{(j>). Modes much greater 
than the characteristic scale qcham are the most signif- 
icantly affected; but as is relevant to the strong cou- 
pling regime we note that as /3 increases, qcham also 



4 



increases. So generally speaking the strongly coupled 
chameleon could evade constraints in yet another way: 
if the scale on which it acts always above observable 
astrophysical scales. 

But it is easy to construct potentials in which this 
is not the case. Recall the fiducial potential has one 
parameter M which has units of mass 



(27) 



For (/) ^ M this reduces to an inverse power law 
V{4>) « M'*+"/(/)". It turns out that the upper bound 
on M that can derived from terrestrial experiments is 
the same energy scale as is relevant to dark energy, 
M < 10"'^ eV 0. A two parameter generalization of 
this, 



= exp( 



(28) 



for an introduction). It is a rather complicated situ- 
ation and much of the astrophysics is poorly under- 
stood and/or constrained. For example is not clear 
whether the satellites are just too faint to be detected 
at present, or if they are actually absent. But surely, a 
persistence of baryons at sub-galactic scales stands to 
worsen the problem, and as Fig. 1 indicates chameleon 
models seem to demonstrate such behavior at strong 
coupling. So it is conceivable that the strong cou- 
pling regime could be constrained by existing astro- 
physical observations, such as the population of dwarf 
spheroidals. This is the topic of future work. 
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has all of the same appealing features as Eq.(27|. If we 
assume M2 ^ (/"min <C Mp;, then 



(29) 



and if we choose Mi ^ IQ^^ eV this is consistent with 
the dark energy scale. It is straightforward to show 
that the slow roll conditions and tracker solution be- 
havior are met as well. While Mi could supply the 
requisite dark energy scale, M2 could operate on scales 
such that it affects the growth of sub-galactic scales 
throughout their evolution. In order to do so we would 
need Pq,j/a > m for roughly the epoch of Zdec to Z},. 
For n — \, these considerations amount to requiring 
13^'^ M2 > 70 eV. 

Such a model could have interesting astrophysical ef- 
fects, for example in the context of the so-called miss- 
ing satellite problem. This phrase refers to the dearth 
of observed satellite galaxies in the vicinity of Milky 
Way sized galaxies compared to what is predicted by 
numerical simulations and analytic estimates (see [13] 
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